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On Small Distuzbances of Plane Couette Flow 1 

Wolfgang Wasow, 2 

The Orr-Sommerfeld equation for the function <f>(y), which appears in the component 
<f>(y)exp\ia(x— ct) \ of the stream function, is analyzed in detail for the case of plane Couette 
flow. A set of solutions is found whose asymptotic behavior as aR— ><» can be calculated 
in the whole complex ?/-plane (R is the Reynolds number). A disturbance with given a. is 
proved stable, for any a, if aR is sufficiently large. In addition, an investigation is made 
of the asymptotic properties of an allied equation for which, in contradistinction to the 
singular Couette case, an inner friction layer does exist. 



1. Introduction 

The purpose of this paper is to analyze mathemat- 
ically some questions connected with the stability of 
plane Couette flow. 

If the space between two infinite parallel planes 
at constant distance in time is filled with a viscous, 
incompressible fluid, and if one plane is at rest while 
the other moves with constant velocity, then the 
simplest flow compatible with the Navier-Stokes 
equations and these boundary conditions is the sta- 
tionary parallel motion with linear velocity profile. 
This is plane Couette flow. Experiments indicate, 
[1, p. 185], 3 that this flow, just as most other viscous 
flows, becomes unstable, that is, turbulent, for high 
Reynolds numbers. 

The oldest and best developed theory of hydro- 
dynamic stability for plane parallel flows, not only 
for Couette flow, studies the stability of disturbances 
so small that they can be regarded as solutions of a 
linearized differential equation. It is further assumed 
that the stream functions of such disturbances can 
be obtained by superposition of functions of the form 



*(y)e ta 



(x-ct) 



(1) 



if the main flow is taken to be parallel to the z-axis. 
Here, a is positive, c and <j>(y) are complex, and t is 
the time. The references [1, 2] contain fuller ac- 
counts of this theory. It reduces the stability prob- 
lem to eigenvalue problems for the differential 
equation 

(4 >-2aV / + «V-iai?[(^-c)((/) ,/ -a 2 ^)-^ ,/ (/>] = O, 

(2) 

where w=w(y) is the velocity of the main flow and R 
denotes a Reynolds number of the flow. 

A solution (eigenf unction) of (2) that satisfies the 
homogeneous boundary conditions of a particular 
flow problem under consideration yields a stable or 
unstable disturbance of that flow according as the 
corresponding eigenvalue c has a negative or positive 
imaginary part. For a given value of R the flow is 

i The preparation of this paper was sponsored (in part) by the Office of Naval 
Research, USN. 
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3 Figures in brackets indicate the literature references at the end of this paper. 



stable, if one has Im c<0 for all possible eigenvalues. 
Otherwise the flow is unstable. 

This approach has been very successful for a 
number of velocity profiles. It has led to critical 
Reynolds numbers that are in good agreement with 
experimental evidence, [4, 5, 6]. A mathematically 
difficult part of this method is the integration of the 
differential equation (2). For large values of aR 
asymptotic methods can be used whose validity 
has been established in [7, 8 and 10]. 

Couette flow, however, although it has the simplest 
velocity profile, plays an exceptional role in these 
theories. First, the proofs of [7, 8 and 10] do not 
cover this case. Second, all investigations seem to 
indicate that Couette flow is stable, in the sense 
specified above, at all Reynolds numbers, although 
no mathematically complete proof of this has been 
given (cf. f2], p. 256). Third, in contrast to other 
profiles, no so-called "inner friction layers" [3, 8] 
occur in the theory of Couette (low. 

This paper will consist of three parts, correspond- 
ing to these three points. 

(a) The first part (section 2) is in the nature of a 
supplement to [9] and [10]. A fundamental system 
of the differential equation (2) for the case of Couette 
flow is constructed whose asymptotic behavior, 
as aR-> oo 9 can be calculated in the whole complex 
y-plane. 

(b) In the second part (section 3) the eigenvalues 
c of (2), corresponding to possible disturbances of 
the flow, are studied as functions of a and aR. It is 
proved that jor any given a a disturbance of type (1) 
of Couette flow is stable, if aR is sufficiently large. 
If the physically plausible conjecture is accepted, 
that above a certain frequency a there is always 
stability (cf. [3], part III, p. 284), the stronger 
statement can be made that, uniformly in a, there 
is a lower bound for aR above which the flow is 
stable. 

On the other hand there are the inequalities of 
Synge ([2], p. 258), which imply the stability of 
Couette flow for sufficiently small values of R. 
Even if taken together, these results do not, of course, 
constitute a complete mathematical proof of stability. 
But the author believes that they go farther towards 
this goal than previous efforts. 

C. Morawetz has shown [13] that the possible 
eigenvalues c of the stability problem fall asymp- 
totically into three classes: (1) Eigenvalues that 
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approach those of the limiting inviscid problem as 
aR-^> oo . The corresponding eigenfunctions approxi- 
mate those of the inviscid problem except possibly 
in certain "friction layers." (2) An infinity of 
eigenvalues lying increasingly close to a certain 
curve in the c-plane as <xR— >oo. They are in no 
way related to the inviscid differential equation. 
The corresponding disturbances are always stable. 
(3) Eigenvalues that approach the points w{y x ) or 
w(y 2 ), as aR->c°, where ij\ and y 2 are the ordinates 
of the flow boundaries. The actual existence or 
stability of eigenvalues of this third kind was not 
discussed in [13]. 

Since no eigenvalues of the first kind exist for 
plane Couette flow, its asymptotic stability, in the 
sense stated above, is proved if the disturbances 
of the third kind are shown to be stable for suffi- 
ciently large values of R. This will be done in 
section 3. 

(c) In all cases examined so far, the actually pos- 
sible disturbances of laminar flows are approxima- 
tions to solutions of the inviscid differential equation, 
except, in general, in certain sectors of the complex 
7/-plane, where they diverge exponentially, as aR-^ oo 
(cf. [8, 9, 10]). The results of section 2 show that 
no such sectors of divergence exist in the Couette 
case. 

In section 4, the mathematical reason for this 
absence is analyzed. Although not immediately 
applicable to physical problems, the results of that 
section may turn out to be of more than purely 
mathematical interest. 

2. Solution of the Differential Equation 

In terms of the dimensionless quantities used in 
(2), the velocity profile of Couette flow is, for an 
appropriate definition of R, 



w(y)=y. 
The further substitutions 

z=y—c, 

\={- a Ri)\ 
<t>(y)=u{z) 
change (2) into the simpler form 



#.W). 



-2a 2 u" + a*U+\ 2 z(u"-a 2 <u) = 0. 



(3) 

(4) 
(5) 
(6) 

(7) 



The theory of [10] does not apply to this differential 
equation, because condition (3.3) of that paper is 
not satisfied. But (7) can be treated directly by 
the method of Laplace transformation used in [9], 
By a straightforward application of Laplace's 
method, it is seen that contour integrals of the form 




FlCURE 1. 

way that &(\"¥)<0. Five such paths are indi- 
cated in the figure 1, where the sectors in which 
Re (\-¥)>0 have been shaded. The paths C(Aj) 
(j = 1,2,3) are the same as in [9 and 10]. The solu- 
tions corresponding to these contours will be called 
Aj(z,X), U(z,X) and V(z,\). Clearly, 



A 1 +A 2 +A 3 =V. 



(9) 



Let S denote any bounded neighborhood of the 
origin in the 2-plane. Then a short calculation 
shows that, in S, 



U(z,\) = cosh az 

a 

V(z,\) = sinh, az . 

a 



(10) 

(ii) 



The asymptotic calculation of the solutions Aj(z,\) 
proceeds exactly as in [9]. It is amply sufficient, 
therefore, to restate those results. 

Let £ and rj be defined by 



HrM-*) 



3/2 



= X 2/3 2, 



(12) 



and denote by C k (&= 1,2,3) the three rays Re £=0 
numbered so that 

4k 2 

arg z=— ir— - arg X, on C k . (13) 



The closed sector of S bounded by Cj^j^k), and 
of central angle 27r/3 may be called S k . If £ is any 
arbitrary positive constant, one has then 



(,,X)-J" o (C-fl-V-""-4>«* (8) (a) Fop ,j|> &| , in s _ ft| 

(M 



are solutions of (7), provided C is a path that is 
either closed or that extends to infinity in such a 



A k (z,X)=^\-™ 



-5/6 



e { [l+0(r)l. (14) 
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Here, the determination of £ for which Re £<0 in S k 
and Re £>0 in S—S k must be taken. 



x-plane 



(b) For|£|<£ , 

A k {z,X) = \- 2/ * ( t-W+'^dt+OQi- 4 *), (15) 

JC*(A k ) 

where C*(A k ) is obtained from C{A k ) by a rotation 
through —2/3 arg X. 

(c) Near G k the asymptotic behavior of A k for 
|{| >; ^o follows immediately from (9), (11) and (14). 

These asymptotic relations may be indefinitely 
formally differentiated. Without much additional 
effort it can be shown that (14) and (15) are first 
approximations of infinite asymptotic series progress- 
ing in powers of £ _1 and X~ 2/3 , respectively. 

The angle of T 5/6 in formula (14) is tt/2- 2tt£/3 
(mod 2tt) at the points of the angular bisector of S— C k . 
At the other points of the domain S—C k , it is then 
obtained by analytic continuation. The angle of 
\-2/3 j s _2/3 arg X. The angle of X will be fixed in 
accordance with (5) by setting 




A convenient fundamental system for the present 
discussion is formed by the functions A u A 2 , U, V. 
Then the eigenvalues are the solutions of the equation 



3w 
arg X=-j- 



(16) F(a,c,\) i 



3. Stability Problem 



M-c) A 2 (-c) U{-c) V(-c) 
Ai(-c) A' 2 {-c) U f {-c) V'(-c) 
Aiil-c) A 2 (l-c) U(l-c) V(l-c) 
A[(l-c)A 2 (l-c) U'{\-c) U'(l-c) 



(22) 



The coordinate system and the definition of the 
Reynolds number can, and will, be chosen so that Two ranges of c have to be distinguished according 
the boundaries of the flow are at as rj c = X 2/ 'c remains bounded or is unbounded as 

X— > °° . Jfrjc is very large, the asymptotic expressions. 
yi = Q, 2/2=1- (17) (14) can be usec [ a t both end points, and arguments 

analogous to those employed in [13] show that the 
The boundary ^conditions for the disturbance are corresponding value of c cannot have positive im- 
aginary part. 

If 7j c remains bounded as X->oo the expressions. 



then ([31, p. 134) 

0(O)=0(1)=0'(O)=^(1) = O. 



(18) 



In terms of the variables introduced in (4), (5) and 
(6) this becomes 

u(—c)=u(l--c)=u'(--c)=u'(l--c)=0. (19) 



(15) are to be used at 2=1- 
(14) for z=— c, that is, r] = r) c 
tion 



-c and the expressions 
Setting, for abbrevia- 



9*W= L,, , t-W^dt (23) 



The differential equation (7) with the boundary one nn d s , after a short calculation, that F(a,c,\)=0 
condition (19) constitutes the eigenvalue problem of becomes in this case 
plane Couette flow. In view of the results of [13] 



described in section 1 only those eigenvalues c=c(a,\) 
need to be studied, for which either 



gi(r, c )+\-™E(a,\,C,Ve)=0 



(24) 



or 



lim c(a,X)=0, 

X->o° 

lim c(a,\) = l. 

X->oo 



(20) 
(21) 



where E(a,\,c,r} c ) is bounded for fixed a, as X— » oo f 
provided 7j c remains bounded. Using Rouche's 
theorem, as in [13], it can be proved that, if 77 = 770 is- 
a root of the equation 



(The notation differs from that of [13] in that there 
X= A /ai?, and that the names of the sectors Si and 
S 2 are interchanged.) The case (20) will be con- 
sidered first. Using (13) and (16), the sectors S k are 
seen to have the appearance indicated by the figure 2. 



then 
that is, 



g'M=o, (25) 

Vc=Vo+0(\- 2 «), (26) 

c(a,X) = -\-^rj + O(\-^). (27) 
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If c(a,X)— »1, the calculations are analogous. The 
roles of the two end points are interchanged, and A 1 
is dominant at the left end point. If 771 is a root of 
the equation 

92(1) =0, (28) 

then there exist solutions of F=0 of the form 

c(a,X) = l-X- 2/3 r 7l + 0(X- 4/3 ). (29) 

The preceding results reduce the question of the 
sign of Im c for large Reynolds numbers to inequali- 
ties concerning the position of the roots of eq (25) 
and (28) . In order to study this question, it is con- 
venient to introduce the variable 



where 

and the function 



S = Q)7), 



iKs)=^-7 gib). 



Since the definition (23) implies that 

.92 (v)=g[ (w»?) = g' 3 (<» 2 v) , 



(30) 
(3D 

(32) 
(33) 



all three functions g' k (rj) can be simply expressed in 
terms of \p(s). 

The values c(a,\) given by (27) have, for suffi- 
ciently large X, positive imaginary parts, if, and 
only if, Re t; >0, that is if the zeros of \//(s) lie in the 
half-plane 7r/6<arg s<77r/6. Similarly, the values 
c(a,\) given by (29) have, for large X positive imagi- 
nary parts, if Re r?i>0. From (33) it follows that 
the necessary and sufficient condition for this to 
happen is that the zeros of \p(s) lie in the half -plane 
— 77r/6<arg s< — 7r/6. Hence all disturbances are 
stable, for sufficiently large Reynolds numbers, if all roots 
of \p{s) lie in the sector |args|<7r/6. 

In order to prove that this is indeed the case, a 
few properties of the function ^(s) must first be 
collected 

From (23), (30), (32) and (33,, it is seen that the 
integral representation 

^( s )= -U gfa)= -U f t-W+Mdt (34) 

YK J 2iri J3K J 27^^ Jc*(A 3 ) 

is valid. By a calculation very similar to that of 
A z (x,\) in [9] asymptotic expressions for \p(s) can be 
derived when |s| is large. The calculation results in 
the following formulas 



*(«) = 



2Vtt 



(-«) 



-3/4 



e 3 



(-S)3/ 2 



[1+0(8-*)] 



(35) 



1_ 1 S -3/*COs(|s 3/2 +^ 



(| S »Hf) 

+ 0(V 7/4 expj| |(Res) 3/2 |}\ (36)" 



Here the first formula is valid outside any sector 
containing the positive real axis, and the second one 
for |arg s|<^27r/3. In the common part of these 
two sectors the two expressions are equivalent. 
In (35) the powers (— s)~* /4 and 2(— s) 3/2 /3 are to be 
taken positive on the negative real axis and con- 
tinued analytically into the whole sector of validity. 
In (36) the powers s~ z/4 and s 3/2 are positive on the 
positive real axis and must be continued analytically 
into |args|<27r/3. The derivation of (35) and (36) 
is omitted except for the remark that actually only 
formula (35) need be proved directly, for from 
formula (33) one obtains analogous formulas for 
g[\rj) and g 2 (y), and the identity 



g' 1 (7l)+g 2 (r 1 )+g' 3 (r J )=27ri, 



(37) 



which follows directly from (23), leads then to (36). 

Differentiation of the relation (34) shows that 

t'(s)=Ai(-s), (38) 

where Ai(s) is Airy's integral (cf. [14] and [8]). By 
(35), the function \p(s) tends to zero as s-^co in the 
sector 2 7r/3<args<4 7r/3. Integrating (38) one finds, 
therefore, 



^( s )= [ Ai(-t)dt. 



(39) 



Where T s is any path starting at infinity in the sector 
27r/3< arg£<47r/3 and ending at t=s. Since the 
Airy integral is real for real s, it follows from (39) 
that \p(s) is also real on the real axis. 

Equations (33) and (37) imply that ^(0) = l/3. 
Hence, (39) may be replaced by 



r 



iKs) = l/3+ Ai(-t)dt 



(40) 



Since Airy's integral is related to Hankel's function 
by the identity 

^ ( ( S )=^Ma(y4 

where the branch of the three- valued Hankel func- 
tion is determined by the condition that Ai[s) must 
be real on the real axis, the integral in (39) can be 
transformed into 



*(«) = 



2iV3 



J L m%(r)dr. 



Here 



Hk 



-sY 



(41) 



(42) 



and i r is a path that begins at r = f and goes to 
infinity in the upper half plane. As the function 
H^Kt) has a branch point at r = 0, while \p(s) is single 
valued, the position of if with respect to the origin 
must be specified. If 



0<args<V, 



(43) 
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and the determination of f in (42) is taken for which 



-7r<argf <-, 



(44) 



a short argument, omitted here, shows that L$ must 
be in — 7r<argr<V/2. 

As to the zeros of ^(s), it follows from the asymp- 
totic representation (35) that outside a sufficiently 
large circle they must all lie arbitrarily close to the 
positive real axis. 

On the other hand, \p(s) has no real zeros. To show 
this, one can use the fact that \p(s) is a solution of 
the differential equation 



^"'+$tf/ = 0. 



(45) 



This can be verified, either by means of (34), or by 
using (39) and the differential equation for Airy's 
function. If (45) is multiplied by 2f, the relation 

d-(2+"+-+' 2 )+s-^m = (46) 

is obtained. Assume, now, that s=a is a zero of \p(s) . 
Integrating (46) from -co to a, it follows then that 



i 



'W- 



tf(t)dt=0. 



This is impossible for real a, because \p(s) is real on 
the real axis. 

After these preparations, it will now be proved that 
the sector 7r/6<arg s< 11/67T contains no zeros of 
\f/(s). Since \[/(s) assumes conjugate values at conju- 
gate points, only the sector 7r/6<args<7r needs to 
be investigated. 

Consider first the portion of this sector for which 



M>r . 



(47) 



If r is chosen so large that the remainder denoted 
by 0(s _1 ) in formula (35) is numerically less than 
unity, for |s| >r , the absence of zeros of \//(s) in the 
domain under consideration is a direct consequence 
of Rouche's theorem. An appraisal of the remainder 
can be based on the following lemma, which is a 
special case of a more general result of Weber [16, 
17, p. 212]. 

Lemma 1. If R(s) is defined by 



then 



H^l(s)=(~y e^~ 5 ^ [l+B(8)] } (48) 



\R(s)\< 



5tt 

72 



11 |- 17 ' 3 / 8 V 

'-iwi ( 1+ 3m) |s|-1 (49) 



for — 7r<arg s<27r. 



Substitution of (48) into (41) followed by an 
integration by parts leads to the expression 

(50) 

From (42), (44) and (47) it follows that f lies in the 
sector — 57r/4<arg f <0 of the r-plane at a distance 
from the origin not less than |r|=2ro /2 /3. The path 
if may be chosen as the ray 



Then 



T={+pe i '< 



<|f" 3/2 ^ f | ^|(r/r) 3/2 ^<-f>|(|+|^(r)r|)|rfr| 
<|r 17 V f l! r - 3/2 (sin|)- 5 / 2 (i+|i?(r*)r*|), (51) 



where 



r*= -^r^sin-- 



For r =9, the expression (51) turns out to be less 
than (0.89)|r iy V r |. 

Reinsertion into (50) shows that the quantity 0(s -1 ) 
in (35) is indeed numerically less than one— even 
less than 0.89 — for 7r/6<arg s<7r,|s| >9. There are 
therefore no zeros of ^(s) in this region. For |«|<r 
the following extension of Rouche's theorem is 
needed. 

Lemma 2. Let \f/(z) be regular analytic in the simply 
connected open region G. Let f(x) be defined and 
piecewise continuous on the rectifiable Jordan curve 
C in G. Assume that on C 



1*00-/(3) |<|/(*) |. 



(52) 



Denote by A arg J the total change of arg f(z) as C is 
once described in positive direction. Then the number 
of zeros of \j/(z) inside C is the integer nearest 
A arg f/2 w. Here arg /(s) can — and must — be defined so 
that it changes continuously at points of continuity of 
f(z) and that it jumps by less than it at the discon- 
tinuities of f(z). 

Proof: The inequality (52) excludes the possi- 
bility of \f/(z) having a zero on C. Let C be defined 
by the equation z=h{t), 0<tf<l. Define a(t) as a 
continuous function such that 

o#)=argiKA(0), 0<KL 

Without loss of generality f(h(t)) can be supposed 
continuous at £=0. Because of (52 J, 



( 1+ ^> 



Re( H 



it. 



273200— 5Sc 
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Hence, 

«(0-arg f(h(f)) 



|arg(l+-^)|<f»'(0<«l) 



for a certain uniquely denned determination of 
arg/(A(£)). Let this determination be called a(t). 
The last inequality implies that a(t) is continuous, 
except for possible jumps of less than w at the dis- 
continuities of f(h(t)), and that 

|[a(0-a(0)]-[a(0-a(0)]K*. 

Since the exact number of roots of ^ inside C is 
[a(t) — a(0)]/2ir, the lemma is proved. 

This lemma will be applied to a contour T in the 
s-plane consisting of three arcs Ti, r 2 , T 3 as follows: 
Ti is the segment from s = to |s|=9 on the ray 
arg s= 71-/6; T 2 is the circular arc denned by s=9e i6 , 
7r/6<0<7r, and T3 is the segment on the real s-axis 
between s=— 9 and 5=0. On r x , the function /(s) 
is to be an appropriate partial sum of the ascending 
power series for ^(s). On T 2 , the leading term of the 
asymptotic series for if/(s) can be taken as f(s). On 
T 3 , the function /(s) is identical with \[/(s). 

The ascending series for \[/(s) follows from (40) 
by term wise integration of the series for Ai(—f), see 
[15] p. B17). Separation of the real and imaginary 
parts for arg s=w/6 produces then two series with 
real terms of alternating sign and decreasing moduli, 
at least from a certain term onward. This fact 
makes it easy to find the number of terms needed 
on this ray for f(s) to satisfy the inequality (52). 

The necessary calculations were first carried out 
at the National Physical Laboratory, Teddington 
(England) under the direction of F. W. Olver. In 
this calculation r was taken equal to 10, and 13 
terms of the series were needed to obtain arg f(s) to 
two decimals of a radian on IV 

Equivalent computations were then performed, by 
a somewhat different method, at the Institute for 
Numerical Analysis of the National Bureau of 
Standards, under the direction of G. Blanch. For 
|s|<6 the values of f(s) on Ti were obtained by 
numerical integration from tables of the Harvard 
Computation Laboratory [18]. For 6<|s|<9 the 
procedure was the same as Olver's. These computa- 
show that the contribution Ai arg / of Ti to the total 
change A arg / of the argument of f(s) is — 3A2ir. 
On his longer ray Olver obtained the value of 
—4.1 It radians. 

On T 2 the inequality (52) is satisfied by virtue of 
the appraisal of the remainder used for | s | >9. The 
contribution A 2 arg / of T 2 to A arg / is 

A a arg /=-§,r+^ 9^=3.43*. 

The corresponding value for |s| = 10 is 4.127T. 

Finally, the contribution A3 arg / is zero, because 
\p(s) is real on the real axis. The condition (52) is 
here satisfied, since \[/(s) was proved to have no zeros 
on the real axis. 



Summarizing, it is seen that A/=0.0l7r on the 
contour formed with the arc |s|=9 as well as on the 
contour with arc |s| = 10 used by Olver. Reference 
to lemma 2 then completes the proof that i/s(s) has 
no zeros in the sector 7r/6<arg s<ll/67r. 

The author expresses his thanks for the advice and 
help he received from Dr. Olver and from Dr. Blanch 
in this section of the paper. 

4. The Differential Equation 
u (4) +X 2 z(" ,r + u)=0 

The occurrence of an inner friction layer is fre- 
quently regarded as a mathematical consequence of 
the fact that the inviscid equation is singular at the 
point where w{y)—c=0. In the case of plane Couette 
flow the inviscid equation, <t>" — a 2 0=O is regular at 
that point, and it seems plausible that this fact is 
responsible for the absence of an inner friction layer 
implied by the unrestricted validity of formulas (10) 
and (11). It will, however, be shown in this section 
that this view is not quite correct. The total absence 
of an inner friction layer is actually a very special 
feature of the differential equation (7). This will 
become apparent from the discussion below, of the 
differential equation 



u™+\ 2 z(u"+u) = 0. 



(53) 



u 
where 



(e,\)= I (<-0*" 1 (i+0" i '" x « ,-+x " a< ' 8/ll " l) *i (54) 
e=l/2i\ 2 . (55) 



Here O must be a contour such that the integral 
exists, and that the contributions from the boundary 
occurring; in the calculations leading to (54) vanish. 
The discussion of these paths in this example is 
made complicated by the fact that the integrand has 



The reduced equation u" -\-u=0 corresponding to 
that differential equation is also regular in the whole 
finite plane. But even a qualitative a priori argu- 
ment shows that some form of an inner friction layer 
is to be expected. The solutions of (53) are all regu- 
lar at z=Q, and their fourth derivatives must there- 
fore vanish there. The reduced equation u"-{-u=0 
has, however, only one linearly independent solution 
with this property. No solution u(z,\) of the full 
equation can, therefore, have a fourth derivative that 
converges to the fourth derivative of, say, cos z in 
a full neighborhood of z=0. But then w(2,X) cannot 
even converge uniformly on any circle | z\ = const, as 
X— >oo ? for u(z,\) is an analytic function of 2 inside 
the circle and it would follow that u (4) (z,\) converges 
uniformly in a neighborhood of e=0, which is a con- 
tradiction. This remark will now be confirmed and 
amplified by actually solving the differential equa- 
tion (53). 

The usual technique of integration by Laplace 
contour integrals leads here to solutions of (53) of 
the form j 
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t -plane 



t -plane 




Ficjure 3. 

branch points at t=±i. In order to bring out the 
analogy with the investigations in [9] seven paths 
are introduced, which are illustrated in the figures 
3 and 4 under the assumption that X>0. 

As to the determination of the multivalued inte- 
grand to be taken the following stipulations are made: 
(1) Two points of the Riemann surface of the inte- 
grand are considered to be on the same sheet, if 
they can be joined by an arc on the Riemann surface 
that does not cross the segment connecting the 
points *=— i and t=i. (2) C(V); C{A k ), (&=1,2,3), 
the lower edge of C(U 3 ) y the upper edge of C(U 2 ), 
and the right edge of C(Ui) lie in the same sheet. 

As before, the functions defined by (54) will be 
called Ate, U k , V according as C is taken to be C(A k ), 
C{U k ) or C(V), respectively. 

From these definitions the following identities can 
be derived: 



A X +A 2 +A 2 =V, 
Uz -U 2 =(e^ 2 -l)A ly 
Ui-Ui=(**-1)A» 



(56) 
(57) 



They result directly from an inspection of the paths 
concerned and the fact that the integrand in (54) 
is multiplied by e r/x2 , if the point t=i is circled once 
in positive direction. The further identity 



Uz-U^^-ViV-A,) 



(59) 



follows by a linear combination of the preceding ones. 
The asymptotic calculation of the A ft (s,X) for 
large X differs only in trivial details from these given 
in [9], and the results are exactly the same. 




Figure 4. 



For V(x,\) one finds without difficulty the asymp- 
totic expression 






(l+* 2 )" 1 «"<ft+0(X-*) 



=2« sin z+0(\- 2 ), 



«;<> 



valid in any bounded domain. 

For the asymptotic calculation of U k for z in S—S k 
and |£|>£ >0, one can proceed essentiallv as in 
[10], p. 12. Since 



x 



t m ci m 

c(u k )l+t 2 dz m 



one has 

few-** 

= f „ T ,[((-i)'" I (*+i)" ( " 1 c x "'W 3 -0-(l+( 2 )- 1 ]re ,2 i/< 



' C(V k ) 



C(U k ) 1 1+t 2 



dt. 



(61) 



The expression in brackets in the last integral is 
numerically less than a constant multiple of 

I X" 2 (f/3 - 1 + arctan, t) | max £ Re [x- 2 (' 3 /3-*+arctan:/] ^ hence 

UC(Uk) 

the integral is of the form 



X ~ 2 Jc(uJ m+l6t ° MX)dt > 



where f(t,\) is uniformly bounded on C(U k ). Split 
C(U k ) into a circle 0\ of radius r about t=i and two 
rays C 2 and C 3 going to, and coming from, infinity 
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at the same angle, in two different sheets of the 
Riemann surface. The contribution of d to the 
integral is 0{\~ 2 ). The contributions of C 2 and (7 3 

k/x \ Hence their sum is 



differ only by the factor — e 



fc 2 



t m + l e tz j(t,z)dt. 



By means of the transformation r=tz, the value of 
this integral is found, similarly as in [10], p. 13, to be 

Oft-V*- 2 ) = 6>(X 2 (w+2) /3 " 4 r 2 (w+2) /3 ) . Hence, 

U£* (z,\) =^ {ire iz ) + OQc 2 ) +0(\ 2 (w+2) /3 " 4 r 2 (m+2) /3 ) , 

(62) 



for z in S—S k and |£ 
in (62; is 0(X~ 2 ) for aJ 



>£o>0,w>0. The remainder 
II m, if s is bounded away from 
zero. Otherwise, the remainder tends to zero for 
m<4, but not necessarily for m^>4. 

In order to deal with the higher derivatives and 
also with the case |£|<£o, observe that the contribu- 
tion to (54) from the path d is in any case equal to 



d m 
dx 7 



(*e iz )+0(\- 2 ). 



The contribution of the two rays C 2 and C 3 is, if a is 
denned by X-¥=<7 3 , 

( e WX 2 _ 1)X 2( M -l)/3 f <e Xp{(T 3 /3- 

X- 4/3 <r+\- 2 arctan a X 2/3 +t,<t}: 



-4/3 



and a short argument, omitted here, shows that this 
is equal to 

X 2(m-l)/3-2 f ^ 3/3+ "V W - 2 rf(T[l+0(X- 2/3 )], 
J C 2 

where C* is a ray in the cr-plane along which Re o- 3 <0 
and whose initial point has at least the distance 
|X~2/3| f rom the origin. 

In S k , finally, the asymptotic behavior of U k (z,\) 
follows immediately from the identities (57) to (59), 
the asymptotic properties of the A k (z,X), and (62). 
It is then seen that XJ k (z,\) diverges for |£| >£ in the 
whole interior of S k in the same fashion as the 
corresponding functions in [9] and [10]* except for a 
factor 6- /x2 -l = 0(X- 2 ). 



Thus, there is in this problem an u inner friction 
layer" of finite thickness for |£| >£ , om^d also a friction 
layer at z=0. The latter occurs only in the derivatives 
of order four or higher. 
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